function outTheta11 = displayResults(theta11,AVarTheta11)

estimates = struct2array(theta11)';
se        = sqrt(diag(AVarTheta11));
tValues   = estimates./se;
% starOne   = abs(tValues) > 1.96;
% starTwo   = abs(tValues) > 2.575;
 starOne   = abs(tValues) > 1.645;
 starTwo   = abs(tValues) > 1.96;
names     = fieldnames(theta11);
for i=1:size(names,1)
    if starOne(i,1) == 1 && starTwo(i,1) == 0
        outTheta11{i,1} = [names{i,1}, ': ', num2str(estimates(i,1)) ,' ',num2str(se(i,1)) ,' * ' ,num2str(tValues(i,1))];
    end
    if starTwo(i,1) == 1 && starTwo(i,1) == 1
        outTheta11{i,1} = [names{i,1}, ': ', num2str(estimates(i,1)) ,' ',num2str(se(i,1)) ,' **' ,num2str(tValues(i,1))];
    end
    if starOne(i,1) == 0 && starTwo(i,1) == 0
        outTheta11{i,1} = [names{i,1}, ': ', num2str(estimates(i,1)) ,' ' ,num2str(se(i,1)) ,'  ' ,num2str(tValues(i,1))];
    end
end


end

